clc
s = 0.7;
hold on
while s < 0.89
    c = [-0.05*(1-s),-0.27*(1-s),-0.19*(1-s),-0.185*(1-s),-0.185*(1-s),s];
    A = [0,0.025,0,0,0,-1
        0,0,0.015,0,0,-1
        0,0,0,0.055,0,-1
        0,0,0,0,0.026,-1];
    b = [0,0,0,0];
    Aeq = [1,1.01,1.02,1.045,1.065,0];
    beq = 1;
    LB = zeros(6,1);
    [x,Q] = linprog(c,A,b,Aeq,beq,LB);
    Q = -Q;
    plot(s,Q,'db');
    s = s + 0.01;
end
xlabel('s'),ylabel('Q')
